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We present the results of numerical calculations on propagation of Extremely-high energy (EHE) 
neutrinos and charged leptons in the earth for trajectories in all phase space of nadir angles. Our 
comprehensive calculation has shown that not only the secondary produced muons but also taus 
survive without decaying in energy range of lOPeV-lOOPeV with intensity approximately three 
orders of magnitude lower than the neutrino flux regardless of EHE neutrino production models. 
They form detectable horizontal or downgoing events in a 1km 3 underground neutrino telescope 
such as the IceCube detector. The event rate and the resultant detectability of EHE signals in 
comparison with the atmospheric muon background are also evaluated. The 90 % C.L. upperlimit 
of EHE neutrino fluxes by a km 2 detection area would be placed at E 2 dF/dE ~ 3.7 x 10~ 8 GeV/cm 2 
sec sr for and 4.6 x 10 -8 for v T with energies of 10 9 GeV in absence of signals with energy-loss 
in a detection volume of lOPeV or greater. 

PACS numbers: 98.70.Sa, 95.85.Ry, 98.70.Vc, 98.80Cq 



I. INTRODUCTION 

It is well known that there exist extremely high en- 
ergy (EHE) particles in the Universe with energies up to 
~ 10 20 eV Q. These EHE cosmic rays (EHECRs) may 
be originated in and/or producing neutrinos by the var- 
ious mechanism. For example collisions of EHECRs and 
CMB photons photoproduce cosmogenic neutrinos 0, 
a consequence from the process known as the Greisen- 
Zatsepin-Kuzmin (GZK) mechanism jSj. Possible pro- 
duction of EHECRs in the present Universe due to the 
annihilation or collapse of topological defects (TDs) such 
as monopoles and/or cosmic strings |^| could also gener- 
ate EHE neutrinos with energies even reaching to the 
GUT scale 0, @- EHE neutrinos provide, therefore, 
an unique probe to explore ultra-high energy Universe, 
which is one of the center piece of high energy neutrino 
astrophysics. 

It has been argued that underground neutrino tele- 
scopes being operated and/or planned to be built are ca- 
pable of detecting such EHE neutrinos • In their travel 
in the earth to the detection volume in a telescope, EHE 
neutrinos collide with nuclei in the rock due to enhance- 
ment of the cross section at EHE range and produce sec- 
ondary leptons like muons and taus. The expected mean 
free path is - 600(/w fc /2.65g cm^-V^/lCT^cm 2 )- 1 
km which is by far shorter than a typical path length of 
the propagation in the earth. Moreover, the decay life- 
time is long enough at EHEs for the produced [i and t to 
survive and possibly reach the detection volume directly. 
Successive reactions of the interactions and decaying are 
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likely to occur in their propagation and the propagation 
processes of EHE particles are rather complex. The accu- 
rate understanding of the EHE neutrino and charged lep- 
ton propagation in the earth is, thus, inevitable for EHE 
neutrino search by underground neutrino telescopes. 

There have been considerable discussions in the lit- 
erature from this point of view. In the Ref. 0, the 
transport equations mainly focusing on v T and r were 
solved and the resultant particle fluxes after the propa- 
gation have been shown for trajectories of several nadir 
angles in the horizontal directions such as 85°. It is true 
that a major faction of EHE r tracks are coming from 
the horizontal directions because the earth is opaque for 
EHE neutrinos, but a km 3 scale neutrino observatory 
like IceCube is essentially a 4-r detector with compara- 
ble sensitivities to both muons and taus, and calcula- 
tion of EHE particle energy spectra of both muons and 
taus over all solid angle space including downward event 
trajectories would be important to evaluate detectabil- 
ity with reasonable accuracy. Furthermore, they utilized 
the often-used continuous energy loss (CEL) approxima- 
tion that follows only the leading cascade particles. It is 
a good approximation for taus, but the secondary par- 
ticle fluxes contributed from the non-leading particles 
are not negligible for muons at EHEs where their de- 
cay does not play a visible role. Calculations on the 
earth-skimming EHE v T have also been made in some 
details They used approximations to neglect contri- 
butions of generated leptons from tau interactions/decay 
in the earth, which would be valid enough for considera- 
tion of earth-skimming neutrino-induced air showers. It 
has been pointed out, however, that the secondary pro- 
duced v e and by tau decays would also enhance the 
total neutrino flux [Toj|. which would be a benefit for an 
underground neutrino observatory. Following all propa- 
gating leptons and taking into account the contributions 
from particles not only skimming but propagating deeper 
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in the earth are, therefore, essential to an underground- 
based neutrino observatory. 

In this work, we numerically calculate the intensity and 
the energy distribution of EHE neutrinos and their sec- 
ondary produced /i's and r's during the propagation in 
the earth for the application to a km 3 scale neutrino ob- 
servatory. The resultant fluxes are shown as a function of 
nadir angles from downward to upward going directions. 
All the relevant interactions are taken into account and 
we follow all produced particles in the reactions whereas 
the CEL approximation follows only the leading cascade 
particles. The initial flux is mainly assumed to be a 
bulk of the cosmogenic neutrinos, generated from the de- 
cay of pions photoproduced by EHE cosmic ray protons 
colliding with the cosmic thermal background photons, 
since the cosmogenic neutrino model is appropriate for 
a benchmark as the flux prediction is on the solid theo- 
retical foundation. Its implications to the detection by 
the IceCube neutrino telescope which are currently 
under construction at Antarctica, are then discussed in 
some details. 

The paper is outlined as follows: First we briefly re- 
view the interactions/decay channels involved with EHE 
particle propagation in the earth in Sec. [H] The method 
of our numerical calculations is also briefly explained. In 
the Sec. IIIII we show the calculated results: energy dis- 
tributions and intensities of muons, taus, and neutrinos 
after their propagation. The energy spectra of these EHE 
particles are shown for the cosmogenic neutrino model. 
Implications on the detection by the IceCube neutrino 
telescope are discussed in the Sec. II VI and the detectabil- 
ity considering the possible background in the experiment 
is discussed in detail. The sensitivity to EHE neutrino 
fluxes by a km 3 neutrino observatory is also shown. We 
summarize our conclusions and make suggestions for fu- 
ture work in Sec. 



II. DYNAMICS OF THE PROPAGATION IN 
THE EARTH 

EHE neutrinos during the propagation do not pen- 
etrate the earth but are involved in charged/neutral 
current interactions that generate charged leptons and 
hadronic showers because their cross sections are ex- 
pected to be enhanced in the ultra-high energy regime. 
Secondary produced /i's and r's travel in the earth initi- 
ating many radiative reactions to lose their energy. The 
higher-order interactions like /i^ pair production [T^| 
and the charged current disappearance reactions like 
fiN — > v^X would regenerate charged leptons and neu- 
trinos which are subject to further interactions. More- 
over, the r's decay channels like r — > v T pLV^ regenerate 
v r . A primary EHE neutrino particle, therefore, results 
in number of particles with various energies and species 
which would be passing through an instrumented volume 
of an underground neutrino telescope. The resultant en- 
ergy spectra and their intensity are consequences from 



TABLE I: Interactions and decay channels involved in the 
EHE particle propagation in the earth. Rows are primary 
and columns are generated particles. 
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the chain processes of interactions and decay. Table |U 
summarizes the interaction/decay channels as a function 
of primary and generated particle species. Main energy 
loss process for secondary produced /i's and r's are, e^- 
pair creation, Bremsstrahlung, and the photonuclear in- 
teractions. The relevant cross sections are formulated in 
Ref. for pair creation, Ref. for Bremsstrahlung, 
Ref. |l5| for photonuclear interaction. Among them the 
photonuclear cross section has the largest theoretical un- 
certainty because it relies on the details of the nuclei 
structure function, which has to be estimated from ex- 
trapolation from the low energy data. In the present 
calculation is used the estimation based on the deep- 
inelastic scattering formalism with the ALLM parame- 
terization of the structure function |l6j | , which has been 
considered to be most reliable prediction. We artificially 
switch off the photonuclear interaction to see its system- 
atic uncertainty in the results later in this paper. Fur- 
thermore, the week interaction, Z iV —* vX, to cause 
muon and tau disappearances, and the heavier lepton 
pair production such as /x + \i~ |l2j are also taken into ac- 
count in the present calculation, which leads to a visible 
contribution to the particle fluxes at EHEs. 

An EHE neutrino is a subject to charged-current (CC) 
and neutral-current (NC) interactions with nucleon. As 
there is no direct measurement of the relevant interac- 
tions in EHE range, the predictions of the vN cross sec- 
tions rely on incompletely tested assumptions about the 
behavior of parton distributions at very small values of 
the momentum fraction x. Since we do not have fur- 
ther clues to investigate EHE neutrino interactions in 
our hands, we limit our present analysis within the range 
of the standard particle physics and use the cross section 
estimated by Ref. 01 using the CTEQ version 5 parton 
distribution functions [l8| . 

Decay processes are also major channels and compete 
with the interaction processes depending on energy. The 
/i and r-leptonic decay distribution can be analytically 
calculated from the decay matrix using the approxima- 
tion that the generated lepton mass is negligible com- 
pared to that of the parent lepton . For z — E Vl / Ei 
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FIG. 1: The energy distribution of EHE leptons after propagation in the earth with nadir angle of 89.5 degree, p's and r's are 
secondary produced. The left panel shows the distributions of leptons with p flavor while the right panel shows the case of r 
flavor. Input spectrum is 10 10 GeV monochromatic of and v T with equal intensity of "1" in this arbitrary unit. 



(/ — p or r) it is written by 



dn 
dz 



3z 2 + -z 3 



3z 2 + -z 3 
3 



(i) 



and for y u = E Ve / (p decay), E Ve ^/E T (t decay), 



dn 
dz 



= 2-6yi + 4y* - (-2 + 12y„ - 18y 2 v + 8y 3 v 



(2) 



The hadronic r decay has various mode and its accu- 
rate treatment is rather difficult. Here we use the 2-body 
decay approximation as in Ref. [2fl |. 

The transport equations to describe the particle prop- 
agation in the earth are given by 



d.J v 
~dX 



= —Na<j v n,cc+nc Jv + 



cprf J Ei dE v 



-:- \ • / dU v ^^£ M e' u ) 
+ .V l/ ,,4^co J , ( ^ ) 



cLh 
dX 



-TT7 = —Na&inJi 



dE v 
mi 
cprfE 



(3) 



Jl 



+Na i dE^^^L J4El) 



+N a I dEi 



dEi 

d(TiN 



dE, 



ME, 



mi 
cprf 



dEi Ir^MEi) 



E[ dEi ' 



(4) 



where Ji — dNi/dEi and J v = dN v /dE v are differen- 
tial fluxes of charged leptons and neutrinos, respectively, 
Na is the Avogadro's number, p is the local density of 
the medium (rock/ice) in the propagation path, a is the 
relevant interaction cross sections, dnf /dE is the energy 
distribution of the decay products which is derived from 
the decay rate per unit energy and given by Eqs.^and|21 
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FIG. 2: Same as Fig.Q but for nadir angle of 70 degree. The 
distributions of leptons with p flavor when the input and 
v T is monochromatic energy of 10 10 GeV are shown. 



c is the speed of light, mi and rf are mass and the de- 
cay life time of the lepton I, respectively. The density 
profile of the rock, p(L), is given by the Preliminary 
Earth Model [2lj . A column density X is defined by 

X = Jp(L)dL'. 
o 

Eq. |21 describes the neutrino propagation. The first 
term is a loss due to the neutrino interaction, the second 
represents a contribution due to the decay, the rest of 
the terms accounts generation of neutrinos by the neu- 
trino and charged lepton interactions. The fourth term 
represents the neutrino appearance by the CC interac- 
tions such as /zN— > i/^X. Eq.0] describes the charged lep- 
ton propagation and has the similar terms with those of 
Eq. 01 but also the term to represent a loss due to the 
lepton decay. 

We numerically calculated these equations by building 
the matrices describing the particle pro pag ation over in- 
finitesimal distance as described in Ref |22ll23l|. The en- 
ergy differential cross sections are derived from the ones 
for the inelasticity parameter y = 1 — E /E, i.e., da/dy. 
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Let us show two examples to show the behavior of the 
EHE particle propagation in the earth. Fig. ^ shows the 
energy distribution of EHE leptons after propagation in 
the earth entering with nadir angle of 89.5°. The cor- 
responding propagation distance in the earth is ~ 110 
km. Primary input spectrum is monochromatic energy 
distribution of 10 10 GeV of and v T with equal inten- 
sity. Sizable bulks of the secondary produced /z's and r's 
are found. As the /x bulk from v T are mainly generated 
from t decay which occurs less frequently in high energy 
region, their intensity decreases with higher energy. For 
the same reason, the secondary r's energy distribution is 
harder than that of //'s. Note that r originated in pri- 
mary denoted as z/ M — > r in the right panel in the figure 
are produced in heavy lepton pair creation fx — > [it + t~ . 

The intensities of "prompt" muon and tau, whose ener- 
gies are approximately same with the primary neutrinos, 
are four to five order of magnitude lower than the primary 
neutrino flux as indicated in the figures, but low energy 
bulk of the secondary muon and tau which has suffered 
energy loss during their propagation makes significant 
contribution to the flux for a given neutrino energy spec- 
trum. It should also be remarked that the muons gen- 
erated from secondary produced tau decay denoted as 
v T — > fi constitute a major fraction of the intensity be- 
low 10 s GeV. We see in the next section that they form a 
sizable flux for the EHE neutrino model producing hard 
energy spectrum like the cosmogenic neutrinos generated 
by the GZK mechanism. 

When particles are propagating more vertically upward 
going, i.e. their propagation distance is longer, all the 
prompt component disappears and no particles essen- 
tially survive in EHE range, because of the significant 
energy losses. A typical case is shown in Fig. [21 for nadir 
angle of 70 degree. One can see that most bulk of the 
secondary muons and neutrinos are absorbed and remain 
only in low energy range. 

Energy distributions and their intensities of EHE par- 
ticles propagating in earth are, consequently, strongly de- 
pendent on the zenith (or nadir) angle of the trajectory, 
and also on the initial neutrino energy spectrum. One 
must solve the transport equation in the entire phase 
space in the zenith angle in order to make accurate es- 
timations of fluxes we see in an underground neutrino 
telescope for a given neutrino initial flux. 



III. THE COSMOGENIC NEUTRINO FLUX AT 
UNDERGROUND DEPTH 

In this section, we discuss the case when the initial 
fluxes of f M and v T are given by the GZK mechanism, 
EHE neutrino production by collisions of EHECRs to 
CMB photons in extragalactic space, as this model has 
been thought to be most conventional mechanism to gen- 
erate EHE neutrinos without new physics and/or spec- 
ulative assumptions. The biggest uncertainty in the in- 
tensity of the cosmogenic fluxes is related to the cosmic 
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FIG. 3: Fluxes of the EHE particles at the IceCube depth for 
a scenario of the neutrino production by the GZK mechanism. 
Two cases in the nadir angle are shown in the figure. 



ray source distributions. Assuming the homogeneously 
distributed astrophysical sources, however, variations of 
the magnitude of the neutrino flux above 10 9 GeV are 
restricted approximately within a factor of 10 |24|. Al- 
though assuming extremely hard cosmic ray injection 
spectrum like ~E _1 or very strong source evolution al- 
lows larger fluxes which can still be consistent with the 
EHECR and the EGRET 7 -ray observations H|, here 
we limit the present calculations to the conventional case 
that the homogeneously distributed astrophysical source 
are responsible for the observed EHECR flux below 10 20 
eV. 

We solve Eqs.|3]and2]to evaluate the particle fluxes at 
an underground depth where a kilometer-scale neutrino 
observatory is expected to be located. The IceCube neu- 
trino telescope is constructed at 1400 m depth below the 
ice surface and we take this number as a representative 
depth. It has been found that changing this depth within 
a factor of two would not affect the overall EHE parti- 
cle intensity in significant manner and the conclusion re- 
mains same. The neutrino oscillation with full mixing 
is assumed and v ^ initial flux is identical to that of v T . 
For t he p arameters constrained by the SuperK experi- 
ment [26| , the oscillation probability in the earth in EHE 
range is negligible, however, and we do not account the 
oscillation in the present calculation on the propagation. 

Fig. |3| shows the fluxes with nadir angle of 85 and 
70 degrees, respectively. The initial primary cosmogenic 
neutrino fluxes are taken from Ref. |2^|. Taus notably 
dominate muons because their heavy mass makes them 
penetrate the earth and because the decay is less impor- 
tant than interactions for the relevant energy range. The 
case of nadir angle of 70° exhibits the strong attenua- 
tion, however, due to the fact that mean free paths of all 
the relevant interactions including the weak interactions 
of neutrinos are by far shorter than the propagation dis- 
tance. This implies that most of the up-going events in 
a neutrino observatory are coming from horizontal direc- 
tions. 

v T flux becomes dominating over that of in low en- 
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FIG. 4: Dependence of the muon and tau fluxes originated 
in the cosmogenic neutrinos on the cosine of zenith angle. 
The integral flux above 10 PeV is plotted in linear scale. The 
atmospheric muon fluxes are also shown by the solid curve for 
the conservative estimation with the low energy extrapolation 
and by the dashed curve for the Corsika-based estimation. 
The detail of the atmospheric fluxes is discussed in Sec. II VI 



ergy range where the r decay is significantly more im- 
portant than interactions. This enhancement is made 
by v T — > t — > v T regeneration process. Note that the 
small bump of the v e spectrum is not the propagation 
effect but generated primarily by EHECR neutron decay 
in the space [HHH]. 

The intensity strongly depends on the nadir angle. 
Fig-EJshows dependences of the secondary muon and tau 
fluxes on the zenith angles. Strong attenuation by the 
earth can be seen but the fluxes are more or less stable 
in the region of the "downward" events where cos 9 > 0. 
Particles in this range are propagating in ice (p = 0.917 
g/cm 3 ) to enter into the detection volume. We numeri- 
cally solved the transport equation in the ice medium to 
derive the downward fluxes. The downward fluxes con- 
stitute major fraction of events in an underground neu- 
trino observatory. The detection issues are discussed in 
Sec. HV| 

The energy spectra integrated over zenith angle are 
shown in Fig. [3] Secondary muons and taus form a po- 
tentially detectable bulk with intensity of ~ three order of 
magnitude lower than the neutrino fluxes. Main energy 
range is 10 PeV to lOEeV (= 10 10 GeV). Regardless of 
the neutrino production model, the relative intensity of \x 
and r to Vy, and v r remains approximately unchanged. It 
should be remarked that intensity of the downward going 
muons and taus are larger than the upward ones by an 
order of magnitude. As also seen in Fig. [31 the tau flux 
dominates over the muons above 10 s GeV. Enhancement 
of v r intensity by the regeneration also appears in the 
upward going trajectories. 

The uncertainty in the muon and tau fluxes estima- 
tions mainly arises from the fact that we do not know 
the photonuclear cross section accurately in EHE range. 
For example, using the updated approach to deduce the 
photonuclear cross section including a soft part of the 



photonuclear interaction would lead to ~ 30 % enhance- 
ment of the total tau energy loss in EHE range p7|. To 
be conservative, in Fig. we show the comparison of the 
fluxes with and without photonuclear reactions. Switch- 
ing off the photonuclear interactions results in a factor 
of two variance in the intensity, which would represent 
the error range of the secondary tau flux estimations in 
a conservative manner. 



IV. DETECTION BY A KILOMETER-SCALE 
NEUTRINO OBSERVATORY 

The event rate for a neutrino observatory can be es- 
timated by integrating of the energy spectra shown in 
Fig. [S] above a threshold energy multiplied by an effec- 
tive area of the detector which is 1 km 2 in case of the Ice- 
Cube. The downward events are major contributions and 
it is necessary to consider the atmospheric muon back- 
ground, however. The atmospheric muon flux estimation 
in the relevant energy range is not straightforward be- 
cause there has been no measurement available and the 
numerical calculation is also time-consuming as one must 
fully simulate EHE air shower cascades. Here we use two 
methods to estimate the flux. One is to extrapolate the 
calculations in 5 TeV which has been confirmed to be 
consistent with the measurement. Because the cosmic ray 
energy spectrum follows E~ 2 ' 7 in TeV region while high 
energy cosmic ray spectrum above 10 PeV are steeper 
following -E" 3 , this extrapolation would be overestimat- 
ing the flux, but it gives the conservative evaluation. An- 
other is to run the Corsika air shower simulation p9l | with 
energy spectrum of the observed E~ 3 under assumption 
that all the mass composition is proton, and count num- 
ber of high energy muons reaching to ground. Then we 
solve the transport equations for the derived muon fluxes 
at surface. The obtained results of the background inten- 
sity in downward events are shown in the right panel of 
Fig. [S] by two dashed lines. One can see that the muon 
background spectrum is quite steep. Setting a higher 
threshold energy, therefore, would be able to eliminate 
the background contamination. The flux dependence on 
the zenith angle is shown in Fig. 0] when the threshold 
energy is 10 PeV. It is clearly seen that the muon back- 
ground attenuates faster than the neutrino-induced EHE 
muons and taus, and there is a window where the signals 
dominate the muon background. Table ^ summarizes 
the intensity with threshold energy of 10 PeV for the 
various EHE neutrino models together with those of the 
atmospheric muon background. 

In fact, what neutrino detectors can measure in direct 
manner is not energy of muon/tau tracks but energy loss 
in a detection volume. The relation between energy and 
energy loss is approximately ~dE/dX ~ f3E. Here (3 is 
the average inelasticity given by 
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FIG. 5: Energy spectra of v e , v^, u T , (jl,t originated in the cosmogenic neutrinos at the IceCube depth. The intensities are 
integrated over solid angle and shown for the upward region (cos# < 0: Left panel) and the downward region (cos 9 > 0: Right 
panel). The two dashed lines represent the atmospheric muon intensities. The upper line shows the conservative estimation 
based on simple extrapolation from the calculation at 5TeV while the lower line is derived by the Monte Carlo simulation with 
the Corsika package. 



TABLE II: Integral flux intensities for several EHE neutrino models. 
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"Intensity at surface before propagating in the earth. 
'Cosmogenic neutrinos with (m,Z max ) = (4.0,4.0) in Ref. \'2'M . 
c Topological Defects scenario using SUSY-based fragmentation 
function in Ref. 0] 
^Estimation based on the Corsika simulation 
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FIG. 6: Dependence of the muon and tau upward fluxes 
on the photonuclear interactions. The integrated fluxes over 
nadir angle of 0° to 90° are shown. 
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FIG. 7: Distribution of energy-loss in propagation of muons 
over 1km in ice. The primary energy of muons is 10 10 GeV. 
Contributions from each interaction are shown separately. 



For the pair creation of muons in ice, /3 ~ 1.3 x 10 6 
cm 2 /g. Therefore the average energy loss fraction due to 
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FIG. 8: Dependence of the muon and tau fluxes originated in 
the cosmogenic neutrinos on the cosine of zenith angle. The 
integral flux above 10 PeV of the energy loss is plotted in 
linear scale. The atmospheric muon fluxes are also shown by 
the solid curve for the conservative estimation with the low 
energy extrapolation and by the dashed curve for the Corsika- 
based estimation. 
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indicating that 10 % of the muon primary energy is de- 
posited in a detection volume. Because the radiative in- 
teractions like Bremsstrahlung have stochastic nature, A 
E is fluctuated significantly in event by event bases, how- 
ever. We carried out a Monte Carlo simulation to see the 
fluctuation. The simulation code uses the same cross 
section and the decay tables but calculates an energy 
of a particle after an infinitesimal propagation length 
AX with the Monte Carlo method instead of solving the 
transport equations. Fig.[7]shows the distribution of the 
energy loss of muons in running over 1 km in ice. The 
energy loss distribution due to the pair creation may be 
narrow enough for the CEL approximation, this is not 
the case for the distributions due to the other interac- 
tion, however. It is not appropriate to approximate the 
entire distribution by a (^-function, which implies that 
the energy loss rather than the energy would be better 
to describe the event characteristics. 

As more realistic criteria, we introduce the threshold 
of the energy loss in ice instead of the energy itself. Fig.0 
shows the GZK integral flux dependences on the zenith 
angle in the case of lOPeV threshold of the energy loss. 
One can see in comparison to Fig.0]that the GZK fluxes 
are larger than or comparable to the muon background 
intensity in all the zenith directions in this energy-loss- 
based criteria. It indicates that it is probable that the 
EHE neutrino search using downward events can be made 
under almost background- free environment. 

It should be noted that the tau flux is lower than the 
muon flux in this criteria. This is because the heavier 



FIG. 9: The tau fluxes at the IceCube depth originated in 
the cosmogenic neutrinos as a function of energy (the dashed 
curve), total energy loss in ice (the solid curve), energy loss 
in form of electromagnetic cascades (the dash-dot curve) , and 
that in form of hadronic cascades (the dotted curve). 



TABLE III: The event rates for several EHE neutrino models. 
The notation of the model name is same as in Table IHl 







,„ > WPeV) 


N T {E U 


,„ > WPeV) 




[km" 2 


year -1 ] 


[km- 2 


year- 1 ] 


GZK Downward 


0.15 




0.10 




GZK Upward 


0.0081 




0.0083 




TD Downward 


1.18 




0.93 




Atm. n 


0.055 








Atm. \i (Corsika) 


0.016 









mass of tau suppresses the energy loss compared to that 
of muons with same energy. This situation is illustrated 
in Fig. where the tau fluxes are plotted as functions of 
energy and energy loss in ice during lkm propagation. 
The intensity above 10 7 GeV is reduced because of the 
energy loss suppression. The higher energy loss takes 
place in form of the hadronic cascades initiated by the 
photonuclear interaction. Table |H] lists the intensity of 
muon and tau above 10 PeV of energy loss for the fluxes 
of the cosmogenic and top-down model @ . The event 
rate under this criteria is found to be 0.27 (fx + r) /km 2 
year for the cosmogenic neutrino fluxes with the moder- 
ate source evolution. Note that the downward event rate 
is 0.25 /km 2 year and dominates in the overall rate. The 
event rates for the various neutrino production models 
are summarized in Table ITTT1 

The IceCube sensitivity on EHE neutrinos can be eval- 
uated by the event rate per energy decade dN/dLogE. 
For a given energy of primary neutrinos, the secondary 
muon and tau fluxes are calculated by the transport equa- 
tions Eqs. and 0] as a function of zenith angles. The 
probability that energy loss with the threshold value or 
greater occurs are estimated by the Monte Carlo simu- 
lation and convoluted with the flux integration over en- 
ergy and zenith angle to give the rate. Fig. ^| shows 
the resultant sensitivity by the IceCube detector with 1 
km 2 detection area. The various model predictions are 
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Energy [GeV] Energy [GeV] 



FIG. 10: The IceCube sensitivities on the EHE neutrino fluxes. 90 % C.L. limits by a 1km 2 detection area with 10 years 
observation are drawn. The left panel shows the case of and the right panel shows the v T case. Labels refer to GZK (|2^| 
for the lower curve, (25| for the upper curve), TD [(J, and Z-burst |3^|. The dashed curves show the sensitivities by events of 
neutrinos interacting inside the detector volume. 



also shown for comparison. The 90 % C.L. upper limit 
i.e., 2.3 event/energy decade/10 year is plotted for 10 
PeV and 1 PeV threshold of the energy loss, respectively. 
The sensitivity is better than that of v T below 10 8 
GeV region because muon energy loss in a detection vol- 
ume is larger than that of taus with same energy, but 
the tau decay which results in large energy deposit in a 
detection volume make dominant contributions in the v T 
sensitivity in this relatively low energy range forming a 
slight bump structure in the sensitivity curve. The 90 
% C.L. upperlimit of EHE neutrino fluxes by a km 2 de- 
tection area would be placed at E 2 dF/dE ~ 3.7 x 10~ 8 
GeV/cm 2 sec sr for and 4.6 x 10~ 8 for v T with ener- 
gies of 10 9 GeV in absence of signals with energy-loss in 
a detection volume of lOPeV or greater. 

This bound would not exclude the cosmogenic neutrino 
production model but strongly constrain the cosmic ray 
injection spectrum in the model. Cosmic ray nucleon 
injection spectra harder than E 1 " 15 would violate the 
bound [2^. On the other hand, as long as the injec- 
tion spectra is softer than E~ 2 , which is very likely in 
case of the astrophysical cosmic ray sources, the IceCube 
bound would less constrain the source evolution as seen 
in Fig. 1101 where we plotted an extreme scenario of the 
cosmological evolution (1+z) 5 where z is the redshift pjj . 
Stronger evolution possibilities than this case are incon- 
sistent with the diffuse background 7-ray observation by 
EGRET [3(| since the GZK mechanism also initiated the 
electromagnetic cascades H3, E3 via photoproduced 
7T° decay and pair creation by EHECRs collisions with 
the CMB photons, forming the photon flux below 100 
GeV which is constrained by the observation. 

The topological defects scenario, on the other hand, 
would be severely constrained by absence of EHE event 
detection by the IceCube. The expected event rate is ~ 2 
events/year km 2 as one can calculate from Table |n] The 
expected EHE neutrino flux in the Z-burst model |3lj . 
the scenario that the collisions of EHE neutrinos with 



the cosmological background neutrinos to explain the 
EHECR fluxes without the GZK cutoff, is well below 
the IceCube bound if the injection neutrino spectrum is 
E~ x as described in Ref. |32l |. 

Although less significant, there are fi and r events pro- 
duced by neutrinos inside the detector instrumented vol- 
ume. In this case the produced charged leptons propa- 
gate only a part of the observation volume. We carried 
out the same Monte Carlo simulation deriving the results 
of Fig. 13 but in which and v r were initially entering 
into the ice volume. The probability that neutrinos in- 
teract inside the 1 km 3 volume and that the produced 
muon or tau losses its energy greater than 10 PeV were 
estimated and convoluted with the neutrino intensity at 
the IceCube depth. The detection sensitivities by this 
channel are shown in thick dashed curves in Fig. 1101 In 
the EHE regime above ~ 10 8 GeV, the intensity of inter- 
nally produced muon and tau events is too small to con- 
tribute the overall sensitivity because the neutrino target 
volume is limited by the size of the detector i.e. 1 km 3 . 
Below 10 8 GeV, on the other hand, including this chan- 
nel improves the sensitivity in sizable manner because 
the energy losses of muons and taus during their prop- 
agation over long distances are more likely to transfer 
them out of the energy range above the 10 PeV thresh- 
old, which leads to reduction of the effective neutrino 
target volume for producing EHE muons and taus out- 
side the detector volume. There is little gain in EHE 
neutrino searches, however, because the proposed EHE 
neutrino models have its main energy range above 10 
GeV. 



V. SUMMARY AND OUTLOOK 

We calculated propagation of the EHE neutrinos and 
charged leptons in the earth to derive their intensities 
and their dependence on nadir angles. The secondary 
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produced muons and talis form detectable fluxes at the 
IceCube depth, with intensity of three order of magni- 
tude lower than the neutrino fluxes. The realistic crite- 
ria, requiring energy deposit grater than 10 PeV in 1 km 3 
volume of ice, leads to ~ 0.27 events/year for the cosmo- 
genic neutrinos in case of the moderate source evolutions. 
The topological defects scenario would be severely con- 
strained. 

The atmospheric muon background are likely to be 
negligible even for downgoing events. The background 
rate is ~ 0.05 event/km 2 year. It should be noted that we 
ignored the possible contributions of prompt muons from 
the charm decay in EHE cosmic ray air showers |33l | , how- 
ever. The atmospheric muon intensity can be increased 
by an order of magnitude but a large uncertainty remains 
due to highly uncertain cross sections in the charm pro- 
duction. The mass composition of cosmic rays in the 
relevant energy range would also be a deciding factor of 
the prompt muon flux intensity: It can be reduced by 
an order of magnitude if the cosmic rays are heavy nu- 
clei |23 and may not constitute a background in the EHE 
neutrino search. Even if the prompt muon intensity is siz- 
able, their energy spectrum would still be much steeper 
than the expected spectrum in the proposed EHE cosmic 
neutrino models, however, and one can easily distinguish 
the signal detections from the prompt muon background 
events if the neutrino observatory has reasonable resolu- 



tion for the energy loss of muons and taus tracks. The 
detector resolution issues require the detailed detector 
Monte Carlo simulations for further investigations. The 
AMANDA experience in relatively low energy muon re- 
construction would lead to energy resolution of A Log E 
~ 0.3 [35j] . The development of the detector Monte Carlo 
simulation is under progress and its application to the 
present results will be an important future work toward 
the search for EHE neutrinos by the IceCube observatory. 
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